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1.   INTRODUCTION 

1.1   GENERAL.   -  -   Marine  navigation  is  the  process 
of  directing  the  movements  of  a  vessel  from  one  point  to 
another.   Historically,  piloting,  dead  reckoning,  and 
celestial  methods  have  been  used  for  this  pi^^oose.   The 
earliest  form  of  marine  navigation  was  piloting,  which 
came  into  being  when  man  guided  his  vessel  with  respect 
to  familiar  landmarks.   Dead  reckoning  o^obably  followed 
as  he  attempted  to  predict  his  future  positions  in  venturing 
beyond  sight  of  land.   While  it  is  generally  accepted  that 
celestial  bodies  were  used  to  steer  by  from  almost  the 
beginning,  celestial  navigation,  as  it  is  known  today,  had 
to  wait  until  information  regarding  the  motion  of  celestial 
bodies  was  acquired. 

The  advent  of  World  War  II  brought  about  the  development 
of  various  navigational  aids  and  electronic  positioning 
systems  which  enabled  the  navigator  to  fix  his  position 
without  regard  to  cloud  cover  or*  weather  conditions.   ^ince 
this  time  the  trend  in  navigation  enuipment  has  been  toward 
obtaining  greater  positional  accuracy,  continous  position- 
fixing  capabilities,  and  longer  range  operations.   These 
improvements  are  being  achieved  by  the  use  of  electromagnetic 
and  acoustic  signals,  coohisticated  electronic  and  celestial 
methods,  and  inertial  systems.   Automatic  computation  has 


become  a  necessity  in  the  more  sophisticated  systems  which 
require  long  and  tedious  calculations. 

The  navigation  requirements  of  the  Navy  are  auite 
different  from  those  of  the  ordinary  user.   An  ideal  system, 
other  than  a  comoletely  self  contained  inert ial  system, 
would  provide  world  wide  coverage,  all  weather  ooeration, 
an  unlimited  number  of  users,  availability  of  frequent 
fixes,  immunity  to  interference,  nassive  user  ooeration, 
and  high  positional  accuracy.   Certain  military  operations 
reouire  oositional  accuracies  of  0.1  nautical  mile  or  better* 

Present  long  range  navigation  systems  and  their 
characteristics  are  presented  in  Table  1.1,    where  long 
range  is  defined  as  having  a  position  fixing  capability  at 
a  range  of  aoproximately  400  nautical  miles  or  greater.   Of 
these,  the  inertial  and  acoustic  doppler  systems  are 
essentially  sophisticated  dead  reckoning  systems  with 
potential  world  coverage  and  position  indicating  capability. 
However,  at  the  oresent  time  and  state  of  the  art,  the  use 
of  a  navigational  satellite  system  appears  as  the  best 
solution  to  the  requirements  outlined  above. 

1.2   CONCEPTION  OF  A  DOPPLER  SATELLITE  NAVIGATION  SYSTEM* 
The  concept  of  a  dopoler  navigational  satellite  system  wa^ 
by-oroduct  of  the  observation  of  the  Soviet  Union's  first 
artificial  earth  satellite,  Sputnik  I,  in  the  autumn  of 


TABLE  1.1  -  LONG-RANGE  NAVIGATION  SYSTEMS 


Estimated  Estimated  Frequency 

System  'TSeful    Range  Accuracy  Range 


Celestial      'Yorld-wide 


3nm 


Star 

Tracker    World-wide 


0.  3nm 


Consol      700-1200nm       0.3-0.6° 


190-1 04k c 
257-263kc 


RDF 


300-600nm 


2-6° 


200-41 5k c 
285-325kc 


lecca 


200-300nm 


0.25-1.0nm     90-1 30k c 


Loran-A    700-900nm 


0.5-3.0nm      1.75-2mc 


Loran-0     1200-1500nm 


G.2-0.5nm     90-110kc 


Omega       5000nm 


0.5-1.0nm  10-1 5k c 


Satellite      World-wide 


0.1-2.0nm     100-400mc 


Ine^tial    3-12  hours 


2.0nm/hr 


Acoustic    300  ft 
Dopoler    water  deoth 


1%   of         0.2-lmc 

velocity 
1%   of  distance 

traveled 


TABLT      1.1       (Cont) 

System  Comments 

Celestial  Inexpensive   equipment,    easy   to   use; 

limited    to    2   fixes/ciay,    affected   by 

clouds  and  bsd  weather  which  also 
limits  fix  availability  to  50°; . 

Star  Tracker        Same  a?  above  except  equipment 

very  expensive. 

Con^ol  Require?  only  a  narrow  band  receiver 

to  establish  LOP. 

RDF  World-wide  coverage.   U.S.  coastal 

beacons.   Capable  of  coarse  fix  only. 

Decca  'Copy  to  use,  insufficient  coverage, 

limited  to  about  240  mile  range  at 
night  because  of  sky  wave  interference. 

Loran-A  Easy  to  use,  insufficient  coverage, 

subject  to  sky  wave  interference. 
Less  accuracy  with  sky  wave. 

Loran-C  Easy  to  use,  accurate,  insufficient 

coverage.   Less  accuracy  with  sky 
wave.   Expensive  installation. 

Ornega  World-wide  coverage.   Presently  in 

advanced  development.   Will  use  three 
frequencies  to  yield  both  eoarse  and 
fine  grids. 

Satellite  World-wide  coverage.   Less  sophisticated 

systems  for  comercial  use  now  unier  study. 

Inertial  Self  contained,  must  be  checked  with 

outside  reference.   Accuracy  decreases 
in  a  few  hours  at  rapid  rate.   Very 
expensive,  requires  expert  maintenance. 

Acoustic  Doppler    Sophisticated  high  accuracy  systems 

utilize  gyro-compass,  velocity  input, 
and  automatic  track-plotting  equipment. 
Limited  to  shallow  water. 


1957.   Dr.  W.  H.  Guier  and  T>.  ft.  q.   Weiffenbach  of  the 
Applied  Physics  Laboratory  of  Johns  Hopkins  Hniversi ty 
listened  to  the  signal  generated  by  Sputnik  and  noted 
the  doppler  shift  in  the  transmitted  radio  frequency.   T'he 
shift  is  given  quantitatively  by: 

(l.i)    Af  =  -(/6/o)f8 

where  A.f  is  the  doppler  shift  from  the  satellite  transmission 
frequency,  y°  is  the  range  rate,  c  is  the  velocity  of 

electromagnetic  propagation,  and  fg  is  the  satellite 

[3] 
transmission  frequency.    Hence,  the  measuring  of  doppler 

shift  is  comparable  to  measuring  the  range  rate.   This  is 

the  basis  for  the  development  of  position  determination 

equations  derived  in  Chapter  3. 

Guier  and  Weiffenbach  deduced,  and  later  demonstrated, 
that  accurate  measurement  of  the  doppler  shift  pattern 
would  permit  the  determination  of  a  satellite  orbit.   After 
the  first  analysis  and  use  of  the  doppler  tracking  technique, 
Dr.  F.  T.  McClure  of  the  Applied  Physics  Laboratory,  reasoned 
that  the  technique  could  be  inverted  and  that  the  doppler 
curve  received  from  a  s  atellite  whose  position  was  known 
could  be  used  to  determine  the  position  of  the  receiver 
on  the  surface  of  the  earth.   It  seemed  only  logical  to 
utilize  this  deduction  as  a  navigation  system. 

1.3   mH^  OPISRATIONAL  RTSQTJIR^M^TT .   -  -   The  T^avy  is 
constantly  seeking  improvements  in  the  fundamental  accuracy 


of  navigation  and  in  the  operational  character  of  navigation 

systems.   Realizing  that  the  doppler  technique  provides  a 

>owerful  tool  for  pimple  'out  very  accurate  navigation,  the 

development  nf  q  highly  accurate  operational  navigation 

system  using  artificial  satellites  was  mode  on  >  of  the  Navy's 

first  sp^C9  objective^.   The  Chief  of  Ha vol  Operations 

established  an  operational  requirement  for  a  navigational 

satellite  as  follow?: 

Develop  a  satellite  system  to  provide 
accurate,  all-weather,  world-wide  navigation 
"ov   naval  surface  ships,  aircraft.,  and 
submarines.  [-4] 

This  reauirement  was  generated  by  the  need  for  a  highly 

accurate  position  determination  required  by  modern  weapons 

systems.   The  satellite  system  was  not  developed  to  replace 

any  of  the  oresent  system?  of  navigation.   Instead,  it  was 

develooed  to  provide  another  navigation  aid  for  determining 

positions  of  extremely  high  accuracy  anywhere  on  the  surface 

of  the  earth. 

Accordingly,  the  Navy  proposed  to  the  Advanced  Research 

Projects  Agency  (ARPA)  that  a  project  be  undertaken  to 

exploit  the  doopler  technique  for  navigational  purposes. 

The  program  began  in  1958  as  a  feasibility  study  under 

contract  with  the  Applied  Physics  Laboratory.   This  program 

was  designated  ''Project  Transit''. 


2.   DEVELOPMENT  OF  AN  OP  NATIONAL  SYSTEM 

•3.1   GENERAL  CONFIGURATION  AND  OPERATION   -  -   After 
the  initial  conception  it  did  not  t'ike  lonp:  to  determine 
the  configuration  of  an  operational  floppier  satellite 
navigation  system.   inhe  basic  system  consists  of  the 
following  five  groups  of  equipment: 

1.  Satellites. 

2.  Tracking  stations. 

3.  Commuting  center. 

4.  Injection  station. 

5.  Navigational  equipment. 

The  system  is  designed  to  employ  four  satellites  in 
pol^r  or  near  polar  orbits  at  an  altitude  of  600  miles  and 
a  period  of  110  minutes.   With  orbital  planes  separated  by 
45  degrees  a  fix  is  available  at  least  every  110  minutes. 
If  the  satellite  oassas  within  10  degrees  of  the  zenith, 
the  interval  is  extended  to  220  minutes  since  longitude 
can  not  be  accurately  determined  in  this  situation.   At 
higher  latitudes  the  freouency  of  fix  availability  increases 
considerably.   Typical  satellite  orbits  for  an  operational 
system  are  illustrated  in  Figure  2.1.    The  satellite 
transmits  two  rigorously  coherent  frequencies  controlled 
from  the  same  oscillator.   This  allows  a  first  order 
correction  to  be  made  for  ionosoheric  refraction.   In 


FIGURE  2.1  -  TYPICAL  OPERATIONAL  SATELLITE  ORBITS 
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addition  the  satellite  transmits  orbital  information  an 
time  signals  synchronized  with  UT-2C.   The  time  signals 
are  obtained  from  a  clock  which  is  regulated  by  the 
satellite  oscillator. 

The  overall  plan  of  a  typical  system  operation  is 
illustrated  schematically  in  Figure  2.2.    t    navigational 
satellite  is  shown  at  three  different  time  periods  during 
an  orbit.   During  time  period  T-,  the  satellite  signals  are 
received  by  a  ground  tracking  station  as  the  satellite 
passes  within  radio  line  of  sight.   As  each  tracking  station 
receives  the  transmitted  signal,  it  determines  the  apparent 
doppler  shift  and  reads  out  the  satellite  memory  and  time 
marker.   These  data  are  then  transmitted  in  digital  form  by 
teletype  to  the  comouter  center. 

After  receiving  a  sufficient  number  of  pass  reports, 
the  comouter  center  uses  the  doopler  data  to  comoute  an 
improver!  satellite  orbit  and  then  Dredicts  the  satellite 
oosition  for  every  two  minutes  in  the  next  16  hour  oerioc. 
When  this  is  completed,  the  computer  generates  the  parameter5 
of  the  orbit  from  which  the  position  of  the  satellite  at 
any  point  in  time  can  be  determined.   The  comouter  also 
analyzes  the  time  errors  to  give  the  clock  rate,  and 
determines  the  correction  to  both  the  clock  setting  and  the 
clock  rate.   All  of  this  information  is  then  sent  to  the 
injection  station. 


During  time  period  Tp,  when  the  satellite  is  within 
range  of  the  injection  station,  the  satellite  memory  is 
erased,  the  new  data  is  inserted  into  the  memory,  and  the 
clock  is  reset  and  regulated.   The  satellite  immediately 
retransmits  this  message  to  the  injection  station  where  it 
is  compared  with  the  information  originally  transmitted.   If 
the  readback  indicates  that  the  injection  was  unsuccessful, 
the  message  is  automatically  retransmitted  to  the  satellite. 
When  accurate  storage  in  the  satellite  memory  is  verified, 
the  memory  is  locked  by  a  time  clock.   Additional  info-:"rr.«ticn 
is  rejected  until  approximately  twelve  hours  later  when 
the  satellite  is  again  in  optimum  range  of  the  injection 
station. 

After  receiving  the  message  from  the  injection  station, 
the  satellite  continuously  transmits  the  parameters  of  its 
orbit  and  its  predicted  orbital  position  at  two  minute 
intervals.   During  time  period  Tg  the  navigation  equipment 
aboard  a  ship  receives  the  orbital  parameters  from  the 
satellite  and  records  the  doonler  shift  for  at  least  three 
2   minute  time  intervals.   Commuting  equipment  uses  this 
information  to  calculate  the  ship's  true  latitude,  longitude 
and  time. 

2. '2   FEASIBILITY  DETERMINATION.   -  -   The  first  phase 
of  the  Transit  Program  was  primarily  concerned  with  proving 
the  feasibility  of  such  a  dopnler  satellite  navigation 
system.   This  phase  was  under  the  direction  of  A^PA. 

11 


During  this  period  four  exDerim<=>ntal  satellites  wore 
designed.   They  were  designated  the  1A,  IP,  ?A,  and  2P 
respectively.   Seven  ground  tracking  stations  were  nut 
into  operation  and  others  were  in  the  process  of  being 
designed  and  constructed.   An  experimental  computation 
center  was  placed  into  operation  and  a  time  and  frenuency 
standards  laboratory  was  established  at  the  Applied  Physics 
Laboratory.   The  requirements  of  the  injection  station 
were  determined  and  the  design  of  the  station  was  undertaken, 
"Extensive  mathematical  analyses  of  computational  procedures, 
ionospheric  refraction,  geodesy,  orbit  determination,  and 
position  determination  were  initiated. 

On  17  September  1959  the  Transit  1A  satellite  was 
launched  by  means  of  a  Thor-Able  vehicle.   The  satellite 
in  it's  flight  shroud  is  shown  in  Figure  2.3.    Although 
the  satellite  failed  to  obtain  orbit,  sufficient  data  were 
obtained  to  verify  the  feasibility  of  s'-telli"Ge  tracking 
and  navigation  by  doppler  analysis.   However,  objectives 
concerned  with  the  verification  of  the  satellite  design 
were  not  achieved. 

The  Transit  IB  satellite  was  similar  to  the  1A  in 
appearance  and  function.   The  satellite  consists  of  a 
shell  divided  into  two  hemispheres  36  inches  in  diameter, 
a  central  supoo^t,  and  an  instrument  tway.   The  shell  is 
a  lamination  of  fiberglas  with  a  honeycomb  plastic  filler. 

12 
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This  provides  a  strong,  highly  insulated,  nonmetallic 

structure.   White  paint  on  the  shell  surface  hida  in 

temperature  control.   The  central  support,  also  made 

of  laminated  fiberglas,  connects  the  hemispheres  of  the 

shell  and  supports  the  aluminum  instrument  tray.   The 

tray  is  tied  to  the  shell  by  nylon  lacing  which  absorbs 

the  shocks  produced  during  launching.   A  cutaway  view 

r  ] 
of  the  Transit  IB  satellite  is  given  in  Figure  2.4. 

The  tray  contains  two  transmitting  systems,  silver- 
zinc  and  nickel-cadmium  batteries,  two  command  receivers, 
and  a  de-spin  system.   A  radiation  and  insulation  shield 
is  placed  above  and  below  the  tray.   Two  banks  of  solar' 
cells  used  for  charging  the  nickel-cadmium  batteries  and 
two  de-spin  weights  are  mounted  outside  and  around  the 
eouator  of  the  shell.   Broadband  antennas  in  a  spiral 
pattern  are  painted  on  the  outside  of  the  shell. 

The  tray  also  contains  an  infrared  scanner  which 
was  designed  for  special  studies  by  the  Naval  Ordnance 
Test  Center".   A  spin  rate  of  approximately  180  rom  was 
imparted  to  the  satellite  for  this  test.   After  a  ore-set 
time  the  de-spin  system  stops  the  spinning  of  the  satellite, 
since  spinning  render's  doppi^r  data  useless. 

The  Transit  IB  satellite  was  successfully  launched 
on  13  April  1960  by  means  of  a  two  stage  Thor- Able -Star 
vehicle.   The  elliptical  orbit  has  an  inclination  of  51.3 
degrees,  a  perigee  of  235  miles,  and  an  apogee  of  475  miles. 

14 
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The  satellite  period  is  95.5  minutes.   Although  the 
satellite  was  not  injected  into  the  planned  circular, 
50O  mile  orbit,  the  ground  tracking  and  computation  oroved 
completely  successful.   Analysis  of  telemetered  data 
indicated  that  the  s at  ell  it e  comnonents  functioned 
Satisfactorily.   All  electrical  and  mechanical  aspects  of 
the  satellite  were  proven  including  the  solar  cell  operation, 
command  off-on  operation,  de-soin  devices,  and  temperature 
stability.   Navigational  experiments  were  conducted  which 
gave  position  determinations  to  within  -1-  mile  of  first  order 
geodetic  control  points.    All  calculations  used  in  orbit 
and  navigational  fix  determination  employed  the  first  order 
ionospheric  refraction  correction  obtained  from  comparing 
the  two  coherent  transmitted  frequencies .   The  accuracy 
obtained  indicated  the  validity  of  this  two  freouency 
technique.   Investigations  of  the  earth's  gravitational  field 
and  geodetic  surface  were  also  conducted. 

The  Transit  IB  satellite  ceased  radiating  on  11  July 
1960.   However,  the  results  of  this  experiment  achieved 
the  goal  of  the  first  phase  of  the  Transit  research  program 
by  conclusively  proving  the  system  feasibility  while 
demonstrating  the  ability  to  perform  dopoler  satellite 
navigation  with  an  accuracy  required  for  most  military 
applications . 

2.5   PROTOTYPE  DEVELOPMENT.   -  -   With  the  success 
of  the  Transit  IB  satellite  the  Transit  Program  entered 
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the  second  ohase  of  development  in  early  Fay  I960  when 
the  management  of  the  Program  was  formally  transferred 
from  ARPA  to  the  Bureau  of  Naval  Weapons.   This  phase 
consisted  of  the  development  of  an  engineering  prototype 
of  the  operational  satellite  navigation  system.   Consequently, 
a  series  of  research  and  development  satellites  were 
constructed  and  placed  into  orbit.   The  basic  characteristics 
and  orbital  information  of  these  satellites  is  given  in 
Table  2.1.    Effort  was  continued  in  improvement  of  system 
design,  simplification,  and  reliability.   Comprehensive 
testing  in  environmental  chambers,  and  for  vibration,  shock, 
and  acceleration  were  conducted  on  the  satellite  and  it's 
components.  -A  complete  data  storage  and  readout  system  was 
developed  including  the  design  and  construction  of  an 
injection  station.   Tracking  orocedures  and  equipment 
were  improved  along  with  the  development  of  shipboard 
navigation  equipment.   Continued  research  was  conducted  in 
the  areas  of  gravity,  refraction,  and  geodesy  resulting  in 
improved  accuracy. 

The  Transit  2A  satellite  was  placed  into  orbit  on 
21  June  1960.   It  differed  from  the  1A  and   IP,  in  that  it  had 
more  solar  cells,  an  improved  telemetering  system,  and  an 
electronic  clock.   The  clock,  running  from  the  satellite 
oscillator,  orovided  the  accurate  time  correlation  between 
the  tracking  stations  which  is  required  for  accurate 
determination  of  satellite  orbits. 
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TABLE  2.1  -  TRANSIT  R^SSARCH  A^:D  DRVKLOPM^NT 
SATELLITES  LAUNCHED  DURING  1960  AND  1961 


*  Note:  Nav  IB  censed  radiating  on  July  11,  1960. 
**  Nav  3B  entered  the  atmosphere  on  March  30,  1961. 
***  Launched  pickaback  on  Transit  4B. 
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The  Transit  3B  satellite,  launched  on  21  February  1961, 
incorporated  an  additional  electronic  unit  in  the  form  of  a 
comolete  data  storage  system  capable  of  storing  a  small 
number  of  bit?  of  digital  information.   A  storage  loading 
technioue  and  readout  system  was  used  to  insure  correct 
data  insertion.   This  provided  the  first  flight  test  of 
the  data  system  including  the  injection  station. 
Unfortunately,  the  3B  satellite  only  remained  in  orbit  a 
little  over  five  weeks  and  the  experimental  results  were 
quite  limited. 

The  Transit  4A  and  4E  satellites  wore  launched  on 
29  June  1961  and  15  November  1961  respectively.   A  cutaway 

re~i 

view  of  the  Transit  4A  satellite  is  given  in  7igure  2.5.- 
The  4  series  was  extensively  engineered  to  achieve  an 
extremely  stable  oscillator,  small  size,  *nd  increased 
reliability.   Crystal  oscillators  of  high  accuracy  were 
carefully  aged  and  selected.   Their  operational  temperature 
was  actively  controlled  by  a  small  heating  coil  contained 
in  the  thermal  inert ial  element.   This  emphasis  on  constant 
operating  temperature  yielded  oscillator  stabilities 
utd  to  one  part  in  10    for*  short  intervals.   Such  accuracy 
gives  the  satellite  system  a  capability  for  extremely 
precise  position  determination. 

The  internal  equipment  was  miniaturized  to  reduce 
the  satellite  size  and  weight.   This  was  to  enable  the 
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use  of  an  inexpensive  Scout  launch  vehicle.   The  outer 
shell  of  solar-  cells  is  collapsable  and  folds  back  igainst 
the  satellite  outer  shell  assembly.   This  greatly  reduces 
the  size  of  the  satellite  as  shown  in  Figure  2.5.   Once 
the  satellite  was  out  into  orbit  the  outer  shell  of  solar 
cells  was  to  be  mechanically  extended.   In  practice  this 
was  never  achieved  since  the  launch  vehicle  was  always  large 
enough  to  allow  the  satellite  to  be  fully  extender!  within 
it's  shroud.   The  4  series  al^o  contained  an  experimental 
radio  isotope  (RIPS)  power  source. 

In  order  to  increase  the  satellite  reliability,  plated 
circuits  and  connection  welding  techniques  were  used  through- 
out.  In  addition,  entire  subsystems  were  imbedded  in  plastic 
to  introduce  the  necessary  ruggedness,  strength,  and 
•temperature  inertia. 

All  the  experimental  satellites  so  far  have  had  to 
use  isotrooic  antennas  because  they  were  not  stabilized 
with  respect  to  the  earth.   Since  no  antenna  is  truly 
isotropic  the  signal  in  some  directions  is  much  weaker 
than  in  others.   As  the  satellite  rotates,  it  sometimes 
sweeps  a  direction  of  low  signal  strength  across  the 
ground  receiver  and  the  signal  lock  can  not  be  maintained. 
This  accounts  for  almost  all  op   the  satellite  pa^sQges 
that  can  not  be  used. 

The  Transit  Research  and  Attitude  Control  (TRAAC) 
satellite  was  launched  pickaback  by    the  same  vehicle  as 
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the  4B  satellite.   It?  primary  objective  was  testing  a 
method  of  gravity  attitude  stabilization.   It  u^es  the 
gradient  of  the  gravitational  field  across  the  dimensions 
of  the  satellite  in  order  tc  keep  one  face  of  the  satellite 
always  toward  che  earth.   This  allows  the  use  of  a  directional 
antenna  and  avoids  the  problem  of  low  signal  strength  in 
the  antenna  pattern. 

As  a  result  of  the  data  gathered  from  the  research 
and  development  satellites  the  basic  characteristics  and 
performance  of  a  satellite  navigation  system  were  validated 
and  a  final  prototype  operational  satellite  was  developed. 
The  first  of  the  operational  Transit  5  series  was  launched 
in  1962.   An  artist's  conception  of  the  Transit  5A  satellite 
is  illustrated  in  Figure  2.6.    Details  pertaining  to  the 
operational  satellite  and  system  are  not  readily  available 
due  to  security  restrictions.   However,  this  satellite 
incorporates  the  design  features  proven  in  the  previous 
experimental  satellites.   Gonspicious  additions  include 
sola7*  cell  paddles  and  a  damping  system.   The  solar  cell 
paddles  can  be  angularly  oriented  with  resoect  to  the  sun 
and  thereby  insure  a  nearly  constant  output  of  power.   The 
damning  system  is  composed  of  a  long  boom  with  a  small 
weight  attached  to  a  spring  at  the  end  of  the  boom.   The 
large  displacements  of  the  spring  provides  the  necessary 
damping  of  angular  oscillations  required  to  keep  the 
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satellite  gravity  stabilized. 

Since  July  1962  the  TT.  S.  Naval  Observatory  has 
maintained  equipment  for  the  recovery  of  time  signals  from 
Transit  type  satellites  to  study  the  problems  associat- 
with  providing  continuous  synchronization  of  UT-2c»    Precise 
time  is  required  for  highly  accurate  orbit  and  position 
determination. 

A  later  model  of  the  -prototype  operational  satellite 
wns  designed  to  use  nuclear  energy  as  a  main  power  source. 
When  this  satellite  was  placed  into  o^bit  it  was  the  first 
operational  satellite  to  use  this  new  energy  in  space.   To 
date,  a  number  of  prototype  satellites  have  been  placed 
into  orbit. 

2.4   THE  OPERATIONAL  SYSTEM.   -  -   In  July  of  1964 
three  of  the  prototype  satellites  were  incorporated  into 
an  operational  system.   A  complete  operational  ground  network 
was  set  up  and  is  entirely  separate  from  the  research  and 
development  network.   The  research  and  development  network 
is  'still  in  existance  and  is  designed  to  provide  basic 
information  on   the  size  and  shape  of  the  earth  and  it's 
gravity  field.   This  information  is  periodically  processed 
through  a  computer  at  the  Applied  Physics  Laboratory  and 
occasionally  new  force  terns  are  derived.   The  new  terms 
are  then  used  by  the  operational  system  to  imorove  orbital 
updating  and  predictions.   In  July  1965  the  name  of  the 
Transit  System  was  officially  changed  to  the  Navy  Navigation 
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Satellite  System.   The  use  of  Transit  is  no  lonper  authorized 
within  the  Navy.   A  ri.ipcupsion  of  the  operational  system 
and  it's  components  is  oresented  in  Chapter  4,  following  an 
investigation  of  the  mathematical  theory  of  satellite 
navigation  presented  in  Chapter  3. 
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3.   NAVIGATIONAL  SATELLITE  THEORY 

3.1   GENERAL.   -  -  We  will  assume  an  artificial 
satellite  in  a  known  polar  or  near  polar  orbit  transmitting 
on  two  extremely  stable  coherent  frequencies  in  trie  100  to 
400  megacycle  range  with  added  uniformly  spaced  time  signals. 
-~  At  a  receiving  station  an  apparent  change  in  the  satellite 
transmission  frequency  is  observed.'  ^his  is  caused  by  the 
motion  of  the  satellite  and  the  receiving  station.   The 
amount  is  prooortional  to  the  velocity  of  approach  or 
recossion.   The  exact  amount  of  this  shift  depends  on  the 
exact  location  of  the  receiving  station  with  respect  to  the 

c 

-J 

path  of  the  satellite.   This  is  illustrated  in  Figure  3.1. 
The  navigation  equipment  integrates  the  heterodyne  frequency, 
which  is  in  the  tens  of  kilocycles  range,  between  the 
satellite  transmission  and  a  reference  oscillator  in  the 
navigation  equipment.   We  note  from  Figure  3.1  that  the 
sign  of  the  frequency  shift  changes  from  plus  to  minus  at 
the  center  of  passage.  '  The  beat  note  does  not  distinguish 
the  sign  of  the  frequency  difference.   In  order  to  get 
correct  measurement^  it  is  required  that  the  beat  frequency 
does  not  pass  through  zero.   Therefore,  the  local  oscillator 
frequency  is  offset  from  the  satellite  oscillator  by  80  parts 
in  10°.   This  offset  exceeds  the  maximum  floppier  =?hift  ^o 
that  the  beat  note  never  changes  sign.   Upon  receipt  of  each 
time  signal  the  number  of  cycles  counted  during  the  previous 

26 


frequency 


time 


OBSERVER  CLOSE 

TO  SUBTRACK 
frequency 

-Af 

yt 

time 


OBSERVER   AT    INTERMEDIATE  DISTANCE   FROM   STTBTRACK 


^  frequency 


+  Af 


-Af 


"Ttime 


OBSERVER  DISTANT  FROM  SUBTRACK 
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time  interval  is  recorded.   The  navigation  eauipment  uses 
the  higher  of  the  two  freauencies  to  determine  count  and  the 
comonri^on  between  the  two  for  a  first  order  correction  for 
ionospheric  refraction.   It  is  then  assumed  that  the  recorded 
cycle  count  is  refraction  free. 

3.2   REFRACTION  CORRECTION.   -  -   Accurate  doppler 
measurements  have  utilized  the  dispersive  nature  of  the 

ionosohere  to  eliminate  the  refraction  contribution  which 

[5] 

is  expressed  as  a  series  expansion  in  the  inverse  frequency. 

f         A(t)    B(t)    C(t)    D(t) 
(3.1)   Af(t,f8)+-i.P(t)=  ~     "  ~<      "T  +  "' 

°  '        fs     fi     fs     fs 

By  receiving  two  coherent  frequencies  transmitted  from  the 
satellite  the  frequencies  may  be  substituted  into  (3.1)  and 
the  first  unknown  coefficient  'A'  determined  at  each  time 
value  t  ,  where  P   is  the  range  rate  and  c  the  velocity  of 
electromagnetic  propagation.   With  the  use  of  only  two 
freauencies  the  resulting  cycle  count  still  contains  all 
contributions  that  depend  on  the  higher  powers  of  (l/f  ). 
In  operational  systems  two  freauencies  *re  generally 
adequate  for  most;  precision  calculations.   The  tracking 
station  receiver  uses  analog  methods  to  automatically  make 
the  refraction  correction  before  the  dopoler  data  are 
digitized.   A  further  di^cus^ion  of  ionosoheric  refraction 
is  oresented  in  Section  5.2. 
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5.5  DETERMINATION  OF  THU  INTEGRATED  GYGL'l  COUNT.  -  - 
The  satellite  transmits  a  time  signal  at  T^  and  at  To.  The 
transmission  time^  from  the  satellite  to  the  receiver  are 

t-|_  -ma    to*   If  we  denote: 

f   -  satellite  transmission  frequency. 
fr  =  reference  frequency. 

f^  =  observed  satellite  transmission  frequency 

'fected  by  doopler  shift. 

IC  —  integrate"!  cycle  count. 
we  can  derive  the  following  relationships. 
The  heterodyne  frequency  observed  at  the  receiver  is: 

(5.2)    th  -  f:  -   f 

N  n     s     r 

and  the  integrated  cycle  count  is  obtained  from: 

Toht 
1   2      2 

(5.5)  IG     =  (fg    -    fr)dt 

^1+tx 

Splitting  the  integral  and  evaluating  the  last  term: 


T2+t2  _ 

(3.4)      IG  -  |       (f*)dt  -  fr  |  (T2+t2)  -   (T1+t1)| 


The  integral  of  (3.4)  is  nothing  more  than  the  number  of 
cycles  transmitted  by  the  satellite  during  the  uniform 
time  interval  Tg  -  Tt. 
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Therefore,    we   can   say: 

m  m 


r 


/• 


(3.5)  I  (fj)dt      =  (fs)dt 

Tlhtl  Tl 

Substituting  the  result  of  (3.5)  into  (5.4): 

(5.6)  IG  =  fs(T2  -  Tx)  -  ^[(Tg+tg)  -  (T^t^l 

By  rearranging  terms  (3.6)  can  be  written: 


(5.7)      IG  =  AfAT  -  f  A" 


Where : 

Af  =  fs  -  fp 

A  T  =  T2  -  T 

At  =  tg  -  t-|_ 

Integration  allows  the  use  of  range  differences  (see 
Equation  3.9)  instead  of  range  rate  in  computations. 

3.4   RELATION  OF  DISTANCE  AND  INTEGRATED  CYCLE  C'TT"T.  • 
The  basic  distance  formula  states  that  distance  is  eoual  to 
velocity  times  time.   Hence,  the  difference  in  transmission 
times  from  the  satellite  is  proportional  to  the  difference 
in  range  at  T-j_  and  T? .   Solving  (3.7)  for  /It: 

Af  AT  -  IG 

(*.B)  At  ='     ~7~ 

1r 
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If  we   ]  et    r    equal    the   range   to    the    satellite  we   car.    express 
(3.3)    in   terms   cf  distance. 

c(AfAT    -   IC) 


(5.9)  A  r    = 


fr 


With  this  basic  equation  relating  distance  and  integrated 

cycle  count  we  can  derive  equations  T"or  position 

determination  in  a  rectangular  coordinate  system. 

3.5   EQUATIONS  FOR  POSITION  DETERMINATION.   -  - 

If  we  denote: 

r  -_  true  distance  to  the  satellite. 

R  -  computed  distance  to  the  satellite  from 
a  DR  position. 

X  ,Y  ,  Z  =  assumed  shio  coordinates  in  a  right 
n   hand  orthogonal  coordinate  system. 

(obtained  from  some  other  navigation 
sv stem) 

Xv,Yv,Zs  —   satellite  coordinates  in  the  same 
l   i   I 

x        system. 

i  =  0,1,2, N  -  where  I!  is  equal  to  the  number 

of  time  oointp. 

Then : 

(3.10)  r  =  R  -f-dR 

Using  assumed  values  for  the  variables  as  opposed  to  the 
true  values  in  (3.9)  and  differentiating  (5.9)  and  (3.10) 
v/ith  respect  to  the  variables  we  arrive  at: 

c 

(3.11)  Ar  =  R1+1    -   R1H-   dR  -    dR_.^   —  (AfnAT    -    ICj 

r 

c  c 

H (ATdAfn)    -—die 
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Expressing  the  R  and  dR  terms  in  terms  of  our  orthogonal 
coordinate  system; 


(3.12) 


R.  - 

1 


OC* 


-  X 


i  -  ):j   (Y'?  -  Y   -  iAy)2 
n       _  x      i    n 


-HZ?  -  7,   -  iAz) 


(3.13) 


dR, 


1 
R. 

-Hz: 


n 


(X? 


n 


-  X   -  iAx)dXr,-j  (Y?  -  Y   -  iAy)dY„ 


-  i;  z)dz 


n 


Where  Ax,  Ay,    and    Az    account    for   navigator  motion  between 
observations.      Similiar    eouations    can   bo   developed    for  R.    , 
and   dR*,,  .      By   substitution  of  values    obtained    above   into 
(3.11)    Ar    can  be    expressed    as: 

(5.14)      Ar  =  Ri+1    -    RjH  j  (X?    -  Xn   -    iAx)dXn   , 


(YT    -   Y      -    iAy)dY        -  (Z?    -   Z      -    iAz)dZ 


n 


n 


n 


n 


-    —       (Xs         -   X      -    (i    l)/.x)dXn -U 

R.     L    1+'1  n 

l-  1 

(Yi+l    •   \   ■    !iJl)i))d'n'1-|ZU    "   Zn   "    U+l)Az)dzJ 

Using   this    expression   for     .'.  r    in    (5.11)    and   rearranging   terms 
we    can    form  an   observation   tvoe    equation  where   v   is    the 
residual    formed    by    the    adjusted    value   of  Ar   minus    the 
observed    value    of   Ar. 


v  =  Arn    -  /\ 
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If  wo  denote  the  coefficients  of  (3.14)  as  follows: 
a.  —  coefficient  of  dx. 


■n 


b^  ■-   coefficient  of  dYn 
c-  =  coefficient  of  dZn 


the  eauation  takes  the  form: 

c 

(.5.15)  v.     —    a.dX     -4-  b    dY    -u   c    dZ AT    dAL-f 

-1  in  i     n  '       i     n        f«  n 


V 


(Vi  -  V  -  -<AfbAT  "  icb' 


In  this  equation  we  assume  our  cycle  count  is  errorless 
and  therefore  neglect  the  dIC  term.   The  navigator  may 
also  assume  that  the  frequency  of  his  reference  oscillator 
is  unknown  but  constant  during  the  short  interval  of 
observation.   In  the  equation  this  is  treated  as  the  unknown 

dAf.   The  above  eauation  contains  four  unknowns  (dX  ,  dY  , 

n»    rr 

dZ  ,  d /if  )  •   If  we  have  four  satellite  observations  we  can 
n'      n 

solve  for  the  unknowns  and  thus  locate  our  position. 

Eauation  (3.15)  has  been  called  an  observation  type 
equation  instead  of  an  observation  equation  because  the 

residual  is  formed  in  terms  of  /\ r  instoad  of  10.   In 
Section  5.7  the  actual  observation  eouation  is  derived 
where  v  is  the  residual  formed  by  the  adjusted  value  of  IG 
minus  the  observed  value  of  IG .   It  will  also  be  shown  that 
for  navigational  purposes  Equation  (3.15)  may  be  used  as  the 
actual  observation  eauation  in  comoutations . 
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If  wo  assume  the  earth  to  be  approximated  by  qn  olliosoid 
of  revolution  then: 


(.6.16)     X2    ■■-  Y;j 
n  '   n 


4- 


__n_ 
b2 


We  can  reduce  the  number  of  unknowns  to  three  by    expressing 

dZ   in  terms  of  dX   and  dY  .   Differentiating  (3.16): 
n  n       n  °  *     ' 


(3.1V) 


d2n=   " 


^n 
—  dX 


Y. 


n 


Z 


n 


-  dY 


n 


1  - 


n 


n 


2 
> 


The  unknowns  are  now  dXn,  dY  ,  and  d/.f.   Now  only  three 
observations  of  a  satellite  are  renuired  to  fix  a  oosition. 
By  substitution  of  the  value  of  dZn  from  (3.17)  into  (3.14), 
multiplying  and  gathering  terms,  we  can  form  an  observation 
equation  similar  to  (3.15). 


(3.18)       V.  =   R      -  R.+ 


1 


Hi  \ 


(Xf  -  Xn  -  iAx) 


:•: 


r       \ 


n. 


(Z|  -  Zn  -  iAzHl-e^X ) 


Zn  .    Ri+1 


(X|+1  -  xn  - 


*_\1 


(ii-l)Ax)-f(Z?4l  -  Zn  -  (i-l)/.z)(l-e2)(-~^)i  dxn 

Zn 


i  (Y? 


v 


R-  * 


n 


-  lAy)+-(zf  - 


y 
Zn  -  iAz)(l-e2)(- 


2n  / 
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1      i 

"—         (YLi-Yn-     (i+D^y)4-(Zj +1-Zn-     (i-KL)Az) 

(i_e2)(__ii)  :  C3Yn (/  e  At  -  ic.)  -  —  At  dz  fn 

Zf»  D  D  « 

n  xr  -J-t. 

This    equation   can  be   reduced    to    the   form: 

(.5.19)  v.    =    a.dX-i-  b.dY    -f    c,d/.f-|-L. 

v  i  in'in'ini 

where  a^  —  coefficient  of  dXn 

bj_  =  coefficient  of  dYn 

c,  =  coefficient  of  d£f 

c   n 

r 
In  practice  we  constrain  our  oosition  to  be  on  a  sohere 

since  we  are  dealing  with  small  Quantities. 

liquation  (3.16)  becomes: 

(3.20)  X^+  Y2-f-  Z2  =  R2 

n    n    n    e 

and  (3.17)  becomes: 

(3.21)  dZn^  -  —  dXn  -  — -  dYn 

Zn        2n 

The  maximum  difference  in  the  dzn  barm  resulting  from 
this  approximation  can  be  determined  as  follows; 

(3.22)  dZn  -  (1  -  &2)dZn 

(ellipsoid)  (sphere) 

and 

dZn  -  ^Zn         ~  -  Q2dZn 

(ellipsoid)      (snhere)  (sphere) 
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Therefore,  the  difference  in  (3.17)  due  to  constraint  to 

o 
'i  sphere  is  -e  'dZn(  pphere)  •   'i'^9  term  dZn  hi?  it's  greatest 

value  at  the  equator.   Assuming  that  the  navigator  position 

is  known  to  on  accuracy  of  one  nautical  mile  (--  1829  meters) 

and  using  the  satellite  derived  value  of  e°  from  the  Kaula 

Ellipsoid  of  1965,  we  find  the  maximum  difference  in  dZn 

to  he; 

AdZn  —  (0.00669)  (1829)   or  ~  12. 2  meters. 

This  result  is  considered  negligible  for  the  purpose  of 
navigation.   By  using  the  spherical  approximation  the 
(1  -  e^)  terms  in  the  coefficients  a^  and  b^  are  eliminated 
and  the  observation  equation  retains  the  form  of  (3.19). 
This  is  the  general  satellite  observation  equation. 

3.6   R"5DTTGTI0TT  OF  OBSERVATION  EQUATIONS.   -  -   Three 
observations  of  a  satellite  will  yield  three  observation 
equations.   In  oractice  these  are  used  directly  in  the 
determination  of  the  three  unknowns.   V/hen  there  are  more 
observations  than  unknowns  the  method  of  least  squares  may 
be  apolied  to  find  that  system  of  values  which  is  most 
probable  from  the  observations  themselves.   Arbitrarily 
assigning  a  weight  of  unity  to  each  observation  and  adding 
the  square  sums  of  the  coefficients  we  obtain: 

(5.23)   [vv]  =  [aa]dx2-t-2  [ab)dXndYn+-  2  [ac]  dXndAfn+2  £>LJ  dXn«h 
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[bb]  dY*  +  2   [bcl  dYndAfn  -  2  |bL]  dYn  +  §c]  dAfg  +  2  [cLj  dAfn-<LL] 

where:   £aa]  =  a^  +  a1+1oi+1   a^ 

By  the  method  of  least  squares  the  most  likely  solution  is 
when  the  sum  of  the  squares  of  the  residuals  (yvj  is  a 
minimum.   To  obtain  this  condition  we  set  the  partial 
derivatives  of  (3.23)  with  respect  to  the  three  variables 
eoual  to  zero.   This  results  in  the  formation  of  three 
normal  equations. 

£a}dXn  +  lab]dYn-hjac]dAfn4-gg  =  0 

(3.24)  Tab|dXn^!bbJ  dYn+jbc]  dAfR4  jbL]  -  0 

|ac]  dXn+  [be]  dYnf  [ccj  dAfn+  \cV   -    0 

"Sauations  (3.24)  are  normal  equations  reduced  from  satellite 
observation  equations.   Regardless  of  the  number  of 
observation  equations  we  can  always  reduce  them  to  a  number 
of  normal  equations  equal  to  the  number  of  unknowns.   The 
solution  of  the  normal  equations  yields  unique  results  for 
the  unknowns  dXn,  dV  ,  and  dAf  •   makinp;  the  values  for  dXn 
and  dY  we  can  solve  for  dZ  .   Adding  these  corrections  to 
the  original  values  for  X  ,  Y  ,  and  Zn  Fives  the  rectangular 
coordinates  of  our  adjusted  position. 

(3.25)  Xa  =  Xn  +  dXn 
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Y  —  Y  4-  dy 

1  a  ~   n  ^   n 

Z  =  Z  -V  dZ 

a     n     n 

The  value  obtained  for  dAfn  may  be  used  to  to  update  the 
freouency  of  the  reference  oscillator  in  the  navigation 
enuipment . 

It  should  be  again  noted  that  the  use  of  the  least 
squares  solution  is  advoided  in  practice  by  using  only 
three  observations  to  solve  directly  for  the  unknowns.   If 
more  than  three  observations  are  obtained  the  excess 
observations  may  be  neglected  or  they  may  be  combined  to 
oroduce  a  total  number  of  three.   This  is  the  most  practical 
solution  ninee  the  advantages  of  the  least  square?  solution 
are  minimized  by  the  small  number  of  observations  that  are 
available  during  one  satellite  pass.   During  a  mcximum  pass 
time  of  approximately  15  minutes  only  seven  observations 
can  be  made.   At  the  present  time  a  decrease  in  the  time 
Interval  to  gain  more  observations  is  not  practical  due  to 
computational  time  requirements  and  restrictions  of  the 
satellite,  memory  and  transmission  system. 

3.7   DERIVATION  OF  THE  OBSERVATION  EQUATION.  -  -   In 
Section  5.5  an  observation  type  equation  was  formed  where 
v  in  the  residual  formed  by  the  adjusted  value  of  Ar  minus 
the  observed  value  of  A^.   In  the  theoretical  case  the 
observation  equation  should  be  derived  in  terms  of  the 
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observed  ouantity.   Such  an  observation  equation  will  now  be 
developed  in  terms  of  the  integrated  cycle  count  where  the 
residual  is  formed  by  the  adjusted  IC  minus  the  observed  IC. 
It  will  then  be  shown  that  for  the  purpose  of  navigation 
the  actual  observation  equation  can  be  replaced  by  Equation 
(3.15) . 
Using  the  notation  of  Section  3.5  and  denoting: 

IC  =  adjusted  integrated  cycle  count. 

IOv.  sss  observed  integrated  cycle  count. 
we  can  express  IC  as  a  function  of  the  unknown  parameters. 


(3.26) 


ICQ  =  f(Xa,  Ya,  Za,  AfJ 


and  : 

(3.27)      IC   =  10.  4-  v    or   v  r~  IC   -  IC, 

v  a      b  '  a      I 

The  adjusted  parameters  can  be  expresses  as: 


v 


Xn4-  * 

yn  +  y 


2a    Zn  |  z 


where  x,  y,  z,  and  F  are  corrections  to  the  assumed  navigator 

parameters.   Assuming  that  the  corrections  are  small  comoarerl 
to  the  parameters  we  can  expand  (3.27)  using  a  Taylor's 
series  where  terms  of  the  second  order  and  greater  are 
neglected . 
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\f  f 

(3.28)      ICh  -t     v    =    f(X    ,    Y    ,    Z    ,  A.f    )    •    -   -x  +  y  -j  z 

D  nnn  n         \  v  vv  '  v 

Xn  -Yn  Zn 


:  f 

r  ....... 

c  Af 

n 

wharv : 

f(x   .   Y  ,   z  ,Af  )   =    IC 

n        n        n  n  n 

The    form   of    the    actual    observation    equation    for   doppier 
satellite   observations    can  be    expressed    asj 

>  f 


(5.29)         v     —     - —  x    ,         -y-t 

^  X  '  Y  r>  Z 

n  n  n 


Z-r   -         F     :-  IC_  IC, 

n 


n 


In  order  to  evaluate  the  coefficients  of  the  corrections  we 
note  from  (3.9)  : 


ic  =  At  A  f  - 


r 


and   by    substitution: 
(3.30)  IC  =  AT£f 


n 


(Xs         -   X    }2  4-    (Ys         -  Y    ^2 
i^i+1        *nJ      ^   ^i+1        *n' 


+  <Zi+l    "Zn> 


(X?   -  xn) 


2 


(Y?    -  Y) 

i  n 


,,i 


<2?    "   Zn) 


The  coefficients  'can  now  be  determined  from  (3.30). 


G 


IC 


^x 


n 


Rw 


(x-Q 
*Ai+i 


I   M-M 


|   /  1 

xj;+  —  (xf 


V] 


OX 


n 
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^1° 


**n 


^Y 


n 


RI+1 


<YLl    "   V 


R, 


<xi   "  V 


r 


dY„ 


_ 


n 


iic 


^z 


n 


Rit-i 


i-hl 


1 


zn) 


fl 


/      iR- 


(Zi    "   V,  J    dZn 


a 


SIC 


=   At  dAf, 
)Af. 


n 


If  we   d  eno  b e  : 


aY  c=  coefficient    of   riX 

bi  — '  coefficient  of  dY 

c*  m  coefficient  of  dZ^ 
1                     n 


we  can  write  (3.29)  in  the  form 


(3.31)         v.     =    a*dXT,  +  b?dy_-V-o*dZn+ATdAf     +   IC_    -   IC^ 

i  i      n        i      n        i      n  n  n  o 

c 

If  we   multiply  '(3.31)    by    a    constant    equal    to we 

f 
arrive    at : 


(3.32) 


c  c 

v,     =   a.rtX_  +  b,dY  4-o,dZ.. ZjT<1  Af 


'i    =   ai'«n+W0i<,Zn 


n 


(icn  -  ICb) 


The  coefficients  of  this  equation  are  equal  to  the 
coefficients  of  (3.15).   From  (3.15)  we  obtain: 

c 


Arn  =  (Ri+]_  -  n±)   =■ 


(ATAfn  -  ICn) 
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and  : 


Arb=     —  (ATAfb  -  ic0) 
v 

(3.15)    can   now   be    expressed    as: 

c 


(3.33)    vi=    aidXn-V-bidYn-i-  CidZn       —  ATdAfn    -       — UCn    -    ICb) 


+  -At  (Af    -  Af.  ) 

f  n  D 

'  r 

From  Equations  (3.32)  and  (3.33)  the  relation  between  the 
actual  observation  equation  and  the  observation  type  equation 
bused  on  terms  of  A r  can  be  expressed  as: 

(3.34)  -  -^-vIG  -     vArV   —  At  (/if     -  Af  ) 

r      i  i  p  n  b 

1r  xr 

and  : 

(3.35)  vTr  =       Kv^°+   KAT    (AfQ    -   Affe) 

c 

where  K  is  a  constant  equal  to  -  . 

f 
r 

Ar  10 

Equation  (3.35)  shows  that  v.    is  related  to  v   by  a 

constant  except  for  the  second  term  on  the  right  side.   It 

will  now  be  shown  that  this  second  term  can  be  neglected 

for  the  purpose  of  navigation.   From  (3.6)  we  note: 
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and  : 


n     s 


rn 


;\ 


f    "   f  V. 

s     rb 


From  this  it  can  be  determined  that  the  term  (Af   -  A fV) 
is  equal  to  (f  ,  -  f   ).   We  will  nov/  determine  the 
aonroximate  maximum  magnitude  of  the  second  term.   Knowing 
that  f  has  a  drift  rate  of  1^««  than  2  parts  in  10    per 
day  and  assuming  that  we  use  the  solution  of  d/. f   from 
previous  observations  to  correct  f   we  can  keen  fr  to 
approximately  2  oarts  in  10-  .   Therefore,  the  term 
( A  f   -  Z  f^)  has  an  approximate  maximum  value  of  (-  2^+-22)i 
or  ±2.8  parts  in  1010.   Using  a  value  of  186,000  miles/sec 
for  c  and  400  mc/sec  for  f   the  approximate  maximum  value 
of  the  second  term  of  (3.35)  can  be  determined. 


KAT(Af  -  A^)  — 
n      b 


186,000 


2.8 


(120)  r-p-  (4  x  10°) 

4  x  10b        101U 


or  dl0.0062  miles  or  £12  meters. 


Since  this  value  can  be  considered  negligible  for 
navigational  purposes  we  can  write: 
Ar 


(3.36) 


Kv 


1C 


Therefore,    if  we    set 


1G     IG ' 
v^    Vj_         to    a  minimum    it    follows    that 
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p      '1 

■  v.'  '  "v.— r   also  goes  to  a  minimum.   Hence  we  are  able  to 

i_  i    i  J 

use  (5.15)  in  place  of  (3.31)  in  cut-  computations. 

5.8   EARTH  RECTANGULAR  CQORDINAT  !S.   -  -   The  adjusted 
coordinate  values  in  a  rectangular  system  are  of  little 
practical  use  to  the  navigatior.   Positional  coordinates 
are  desired  in  terms  of  latitude  and  longitude.   While 
the  calculations  for  coordinate  determination  and 
transformation  are  generally  well  known,  they  are  included 
to  provide  a  complete  development  of  a  position  determination 
from  donpler  satellite  observations. 

In  calculating  the  positions  of  navigational  fixes  on 
the  surface  of  the  earth  it  is  necessary  to  consider  these 
points  as  lying  uoon  some  mathematical  surface.   The  actual 
shape  of  the  earth's  surface  is  quite  irregular  and  is  not 
adapted  to  the  purpose  of  computation.   For  this  reason  it 
is  necessary  to  adopt  an  assumed  figure  which  will  not  d  ^Dart 
from  the  true  surface  by  an  amount  sufficient  to  produce 
serious  errors  in  our  results.   The  figure  generally  adopted 
is  an  ellipsoid  of  revolution  which  is  defined  by  a  given 
flattening  (f)  and  a  given  semi-major  axis  (a). 

The  position  derivation  given  in  the  preceeding  sections 
is  with  respect  to  an  arbitrary  right  hand  orthogonal 
coordinate  system.   We  will  now  define  this  system  with 
resoect  to  the  ellipsoid  of  revolution  assumed  to  represent 
the  figure  of  the  earth.   The  origin  is  at  the  center  of  the 
ellipsoid  which  is  also  the  assumed  center  of  mass.   The 

45 


Z  axis  coincides  with  the  mean  rotational  axis  an'']  the  X 
axis  intersects  the  Greenwich  Meridian  at  the  equator.   The 
Y  axis  completes  the  right  hand  system.   The  coordinate 
system  thus  defined  is  illustrated  in  Figure  5.5. 

The  relations  between  geodetic  latitude  and  longitude 
and  the  defined  XYZ  coordinate  system,  which  can  be  derived 
from  Figure  5.5,  are  given  in  the  following  formulas. 


X  = 


(5.57)        Y 


a  cos  (p  cos  A 

W 
a  cos  Cp  sin  (\   , 

W 
a(l  -  e2)  sin(p 


'Vhore  W  is  defined  'is: 

W  =  (1  -  Qd    sin"  Cp  ) 


5.9   DETERMINATION  OF  LAmITT^)v  &ND  LONGITUDE.  -  - 

Using  values  Tor  X  ,  Y  ,  and  Z   from  (5.25)  and  the  relations 

c  a    a        a 

given  in  (5.57)  we  are  able  to  determine  the  latitude  and 
longitude  on  the  surface  of  the  given  reference  ellipsoid. 
The  value  for  longitude  can  be  solved  directly  from  (5.57). 


(5.58)      A  =  Tan-1  

X 
a 
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n  Section 


a 
b 

e 


semi  major  axis 
semi  minor  axis 
first  eccentricity 
geodetic  latitude 
geodetic  longitude 


4T. 


a 


FIGURE  3.3  -  bARTH  RECTANGTTLAR  COORDINATE  SYSTEM 
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The  solution  for  latitude  is  s  little  more  involved.   From 

Figure  5.4  we  observe  that: 

N  sinCp 
(5.59)     tan(p  =        x 

(X2  Y^s 

Where  N  is  the  radius  of  curvature  of  the  surface  of  the 
ellipsoid  in  a  plane  through  the  normal  and  at  right  angles 
to  the  meridian.   This  is  also  referred  to  as  the  radius 
of  curvature  in  the  prime  vertical.  A   problem  arises  that 
N  is  a  function  of  Q)  .   (N  is  equal  to  a/W)   To  solve 
(5.59)  we  must  eliminate  (0   from  the  right  side  of  the 
equation.   Going  back  to  (5.57)  we  can  express  Z  in  the 
following  manner: 

(5.40)  Z  =  N  sin  (p  -  N  e2  sin  Q) 

and  N  sin  0  —    Z  -\-  N  e2  sinQp 

Substituting  (5.40)  into  (5.59): 

ZfN  e2  sinCp       Z        N  e2  sinCp 

(5.41)  tan(T)~  — — =         1  (1+ ) 

(X2  Y2)2        (X*-  Y2)3        Z 

Substituting  N(l  -  e2)  sinQj)   for  Z  in  the  denominator 
of  (5.41)  we  arrive  at: 

-1      Z* 


(5.42)      (p  =  Tan""L( —  ) 


(X2  Y2)^   1  -  e2 
a   a 

The  values  of  latitude  and  longitude  obtained  from  (5.58) 
and  (5.42)  are  used  by  the  navigator  to  fix  his  position 
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p  1  P" 

(X+Y^)a 


FIGURE  3.4  -  MERIDIAN  SECTION  IN  EARTH  COORDINATE 
SYSTEM. 
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on  the  surface  of  the  ellipsoid. 

Reversing  this  procedure  the  navigator  can  easily 
calculate  his  position  in  earth  rectangular  coordinates 
from  values  of  latitude  and  longitude  by  applying  (3.37). 

3.10  INERTIAL  COORDINATE  SYSTEM.   -  -   The  derivations 
to  this  point  have  assumed  that  the  navigator  is  given 
satellite  coordinates  in  the  defined  earth  rectangular 
coordinate  system.   In  practice  the  navigator  must  calculate 
these  coordinates  from  given  orbital  parameters  which  give 
satellite  coordinates  in  an  initial  system.   The 
determination  of  these  parameters  is  quite  complex  and  is 
omitted  from  this  thesis.   Renders  desiring  this  information 
are  referred  to  Reference  7. 

The  defined  inertial  coordinate  system  is  illustrated 
in  Figure  5.5.   The  origin  of  this  system  is  the  same  as 
that  of  the  earth  rectangular  coordinate  system  defined  in 
Section  3.8.   The  Zj  axis  coincides  with  the  pole  of  the 
orbital  plane.   The  Xj  axis  passes  through  the  perigee 
and  the  YT  axis  completes  the  right  hand  system. 

3.11  SATELLITE  COORDINATE  DETERMINATION.   -  -   All 
navigational  satellite  systems  require  accurate  orbital 
information  so  that  the  satellite  coordinates  at  any  point 
in  time  may  be  calculated.   Since  it  is  not  possible  to 
predict  the  position  of  a  satellite  with  respect  to  time 

to  the  required  accuracy  for  more  than  a  day  it  is  impo^^ible 
to  include  a  satellite  ephemeris  in  a  nautical  almanac. 
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I  H 


FIGURE  3.5  -  INERTIAL  COORT)INAT^]  SYSTEM 
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To  overcome  this  drawback  tho  satellite  contains  an  ephemeris 
table  stored  in  a  magnetic  memory  which  is  periodically 
updated  with  new  parameters  from  the  ground. 

The  memory  is  divided  into  two  parts.   The  first  or 
fixed  Dart  contains  the  oarameters  required  to  define  the 
best  elliotical  approximation  to  the  actual  orbit.   These 
oarameters  change  slowly  with  time  and  are  taken  as  constant 
over  a  twelve  to  sixteen  hour  period.   The  discrepancies 
between  this  orbit  and  the  actual  orbit  are  small.   The 
second  or  ephemeral  part  of  the  memory  contains  the  small 
corrections  to  the  basic  orbit  for  each  time  point.   The 
orbital  parameters  stored  by  the  satellite  memory  are  given 
in  Table  3.1.   From  this  data  the  navigator  is  able  to  compute 
the  position  of  the  satellite  in  inertial  rectangular 
coordinates  at  the  beginning  and  end  of  each  time  interval. 

To  calculate  the  satellite  coordinates  x!?,  Y?,  and  Zs 
the  navigator  goes  sequentially  through  the  following 
calculations  where  tm  is  the  mean  satellite  pass  time. 

(3.43)       tm  =  tm  -  tp 


(o),  /  \     At, 


pP    =  (pPv°VuJp- 

AN  =  Aa(o)+W*At» 


Mm  =  M0D(h.AM  +  GJo£tm 

ooO 
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TABL^   3.1  -  SATELLITE  ORBITAL  PARAMET^S 


Fixed  parameters: 

t    a  time  of  first  perigee  after  0h  or  12hTTT. 

»     v  mean  angular  motion  «  2  "a/period. 

(o) 

#  argument  of  perigee  at  t  . 

|jj    *  precession  rate  of  perigee. 

C_  -^  eccentricity  of  orbit. 


AQ   *  mean  semi -major  axis, 

A  (°) 

_y\_AT   *  right  ascension  of  ascending  node  at  t 
K  P 


Ll«  *  precession  rate  of  ascending  node. 

C;  J«-  cos  {   *  ^inclination  of  orbital  plane 

S(  *  sinj. 

Am  *  change  of  mean  anomaly  for  1  hour. 

<£  M  *  change  of  mean  anomaly  for  2  minutes, 
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TABLE   3.1  (Cont) 

Ephemeral  parameters: 

t^   *   time  nfter  intepral  half  hour   ^.,  A*. 

£ E •  *   correction  to  true  anomaly  at  t^. 

<3  A.  *   correction  bo  semi-major  axis  at  t^. 
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M 


i  =  *m  +  i\  -  fcJ  Sm/. 


6*i    =    £>in  Mi 

Ai    -    Ao^Ai 

Then  the  satellite  coordinates  in  the  inertial  system  at 
time  t^  -ire  fiven  by: 

X*  =  Aj_(  cos  E±    -  £_) 

(3.44)      Y^  =  a±  [sin  7^(1  -  cE-2)i] 

Zj  —  0 

Theso  are  the  satellite  coordinates  in  the  defined  inertial 
coordinate  system  at  the  i  th  time  point.   In  order  to  be 
used  for  position  determination  in  (3.18)  the  inertial  co- 
ordinates must  be  transformed  to  the  defined  earth 
rectangular  coordinate  system. 

3.12   TRANSFORMATION  OF  INERTIAL  COORDINATES.   -  -  The 
transformation  of  the  inertial  satellite  coordinates  to  the 
earth  rectangular  coordinate  system  can  be  considered  in 
two  distinct  steps.   The  first  step  involves  two  rotations 
to  bring  the  inertial  XT  and  Yy  axes  into  the  equatorial 
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plane  and  the  Zy  axis  to  coincide  with  the  rotational  axis. 

t    t        . 

The  results  of  these  rotations  define  the  XTf  Y  ,  and  Z 

I'  1'  i 

axes  illustrated  in  Figure  '5.5,      The  second  step  involves 
bringing  the  X-r  axis  to  coincide  with  the  X  axis.   This  is 
complicated  by  the  earth's  rotation  about  its  axis*   Since 
the  earth  rectangular  coordinate  system  rotates  with  respect 
to  the  inertial  coordinate  system  the  second  step  takes  the 
form  of  a  rotation  of  the  inertial  prime  system  about  the 
Z  axis  through  an  angle  of  - -A-w  ~  GAST .   Upon  completion 
of  these  steps  the  satellite  coordinates  will  be  in  the 
earth  rectangular  coordinate  system  used  in  the  position 
derivation.   The  matrix  form  of  the  required  rotations  is 
given  below; 


Step  1  -  rotation  about  the  Zy  axis  through  an  angle 

„(o) 

P 


of  -  y     followed  by-  a  rotation  about  the 


Xj  axis  through  an  angle  of  -  ( 


(3.31) 


r  r 


z! 


0 


0 


0   cos(-()   sin(-0 


0  -sin( -I)      cos(-(  ) 


cos  (  - 


(o) 


,.Ao) 
-sin(-<y 


0 


(o)    -^  ,- 

sin(-  (P   )   0  !  X. 

p        I 


cos(-Cpp   )   0 


0 


y. 


l]  j  z 
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Step  2   -   rotation  about  the  Z  axis  through  an  angle 


of  (-  JVN   -    GAST; 


(3.32) 


i_  _t 


cos(AN    -    GAST) 


sin(AN   -   GAST) 


0 


•sin(AN   -   GAST) 


cos(AN   "    G^ST) 


0 


0 

1 


1 

yI 


LziJ 


3.13   SUMMARY   -  -   The  theoretical  development  of 

eouation^  for  position  determi nation  given  in  the  prec'eeding 

sections  illustrates  how  a  navigator  is  able  to  fix  his 

position  on  the  surface  of  the  ellippoid  user)  as  nn  assumed 

figure  of  the  earth  from  observations  of  a  navigational 

satellite.   Initial  assumed  coordinates  X« ,  "V    and  % 

n 7   n       n 

are  reonired  from  some  independent  source.   In  addition, 
the  navigator  must  determine  the  integrated  cycle  count 
from  at  least  three  satellite  observations  and  be  able  to 
calculate  the  coordinates  of  the  satellite  at  the  time  of 
these  observations.   The  solution  involves  the  determination 
of  corrections  dX  ,  dYn,  and  dZ   which  are  applied  to  the 
assumed  initial  coordinates.   This  gives  the  new  position 
coordinates  in  a  rectangular  system.   These  coordinates 
are  then  expressed  in  term^  of  geodetic  latitude  and 
longitude . 
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4.   OPERATIONAL  SYSTEM  COMPONENTS 

4.1   GENERAL.   -  -   The  Navy  Navigation  Satellite 
System  has  been  operational  since  July  1964.   It  consists 
of  three  satellites,  four  tracking  stations,  an  operations/ 
computer  center,  two  injection  facilities,  and  numerous 
navigational  receiving  units.   The  ground  network  is  based 
entirely  within  the  United  States  and  is  operated  by  the 
U.  S.  Navy  Astronautics  Group.   The  locations  of  the 

M 

components  of  the  ground  network  are  shown  in  Figure  4.1, 

The  task  of  the  Astronautics  Group  is  to  determine, 
process,  and  load  orbital  and  time  information  into  the 
memories  of  the  operational  satellites.   These  operations 
must  be  carried  out  within  critical  time  limits.   Each 
satellite  is  tracked  every  time  it  comes  within  range  of 
one  of  the  four  tracking  stations.   The  data  received 
from  the  satellite  is  then  sent  to  the  computer  which 
updates  the  orbit,  predicts  the  satellite  positions  for 
the  next  16  hours,  and  generates  the  orbital  parameters 
in  a  coded  message.   This  message  must  be  sent  to  the 
injection  facility  in  sufficient  time  to  be  injected  into 
the  satellite  before  it  has  used  up  all  of  the  information 
previously  stored  in  its  memory.   Each  cycle  requires 
about  12  hours  per  satellite.   Since  the  satellite  memory 
is  limited  to  16  hours  of  information,  operations  must 
meet  an  exacting  schedule. 
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4.2  THE  SATELLITE.   -  -   As  was  previously  mentioned, 
details  concerning  the  operational  satellite  are  not  readily 
available  due  to  security  restrictions.   However,  it  can  be 
assumed  that  the  operational  satellite  employs  the  basic 
characteristics  developed  by  the  Transit  research  and 
development  series  described  in  Chanter  2. 

4.3  TITS  TRACKING  STATION.   -  -   The  four  operational 
tracking  stations  are  located  in  lAaine,  Minnesota,  California, 
and  Hawaii.   The  tracking  station  equipment  consists  of  a 
Westinghouse  AN/BRN-3  type  receiver,  a  doppler  cycle  counter, 
an  ultra-stable  oscillator,  and  communications  equipment. 

The  data  obtained  by  the  tracking  station  is  put  into  standard 
teletype  format  consisting  of  two  parts.   The  first  part  is 
called  the  header  and  contains  the  identifying  information 
required  by  the  computer.   The  second  part  contains  the 
actual  data.   The  data  consists  of  a  series  of  11 -digit 
data  points.   ^ach  ooint  represent?  one  doppler  reading  and 
the  time  that  the  reading  was  taken.   These  readings  are 
usually  taken  every  four  seconds.   Generally  100  or  more 
data  points  are  obtained  during  the  pass  of  one  satellite 
within  radio  line-of-sight  of  the  tracking  station.   In 
messages  sent  to  the  operations/computer  center,  the  first 
five  digits  of  a  data  point  represent  the  time  of  day  in 
seconds  (  U  T  )  and  the  last  six  digits  represent  the  number 
of  microseconds  it  takes  to  receive  a  predetermined  number 
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of  cycles  from  the  satellite.   This  preselected  value  is 
such  that  the  time  approaches  but  never  equals  one  second. 

In  order  to  send  the  great  quantity  of  data  gathered 
by  the  tracking  stations  to  the  ooerations/computer  center, 
high  speed  communication  eauioment  is  required.   A  rapid 
data-link  called  TRADAT,  which  has  a  transmission  rate 
five  times  faster  than  teletype,  is  presently  being  used. 
TRADAT  is  essentially  error  free  because  it  incorporates 
self-checking  circuitry. 

4.4   OPERATIONS/COMPUTER  CENTER.   -  -   The  operations/ 
computer  center  is  located  at  the  headquarters  of  the 
U.  5>.  Navy  Astronautics  Group  at  Point  Mugu,  California. 
All  data  from  the  tracking  stations  and  the  computer  pass 
through  the  operations  center.   This  center  at  all  times 
maintains  the  operational  system  status  and  the  positions 
of  the  satellites  as  well  as  related  systems  information. 
It  is  the  operations  center's  responsibility  to  insure 
that  the  tracking  stations  track  each  satellite  and  that 
sufficient  data  is  provided  to  the  computer  to  allow 
accurate  orbit  determ1  lations  and  predictions.   It  schedules 
the  computer  center  runs  and  insures  that  the  orbital  and 
time  parameters  are  sent  to  the  injection  facility  in  time 
to  meet  the  injection  pass  of  a  satellite. 

Four  satellites  in  orbit  will  produce  an  average  of 

84  messages  per  day  from  the  tracking  stations.   The  data 

in  these  messages  pass  through  the  operations  center  and 
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are  sent  bo  bhe  computer  center  whose  primary  ^oo    is 
orbit  updating  and  position  prediction.  The    computer 
center  has  an  IBM  7094  II/1410  configuration  with  a 
360/40  for  command  and  control  functions.   The  dcopler 
data  are  fed  into  the  computer  which  is  programmed  for 
orbital  analysis.   The  machine  computes  a  theoretical 
orbit  of  the  satellite  and  compares  a  theoretical  doppler 
curve  obtained  from  this  orbit  v/ith  the  actual  doopler 
curve.   The  residuals  are  squared  and  summed  and  a  new 
set  of  corrected  variables  is  obtained.   These  allow  the 
computation  of  a  new  theoretical  orbit  and  a  new  theoretical! 
doppler  carve  to  be  again  compared  v/ith  the  observed  data. 
This  process  is  repeated  until  the  sum  of  the  squares  of 
the  residuals  is  a  minimum.   .At  this  point  it  is  assumed 
that  the  best  determination  of  the  orbit  has  been  made. 
In  this  way  orbits  of  high  precision  can  be  determined 
and  the  future  position  of  the  satellite  car.  be  predicte 
with  bhe  required  accuracy.   It  takes  nearly  two  hours 
to  update  a  single  satellite  orbit  and  prepare  an  injection 
message.   V/ith  each  satellite  requiring  this  information 
approximately  every  twelve  hours,  the  computer  will  be  used 

to  capacity.   The  message  prepared  by  the  computer  contains 
:4,960  bits  of  orbital  and  time  information  which  comprise 

the  message  proper  and  320,000  bits  of  information  which 
control  the  equipment  at  the  injection  facility  during  a 
satellite  pass. 
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^Ithouph  TRADAT  speed  and  accuracy  are  particularly 
i&apted  to  communications  between  the  tracking  stations 
>:n.d    the  operations/computer  center,  it  still  does  not 
provide  the  speed  and  accuracy  level  required  to  handle 
the  injection  messages  which  are  sent  from  the  operations/ 
computer  cent<  ?   to  the  injection  facility.   These  must  be 
100  per  cent  error  Pre<  .   High-speed  equipment  with  a  data 
rate  more  than  twice  that  of  TRADAT  is  used  for  the  sole 
purpose  of  sending  injection  messages,  which,  from  the 
time  they  are  fed  into  the  equipment,  are  entirely  automated 
until  they  are  transmitted  by  the  injection  facility 
a nt enna  to  the  satellite. 

4.5   THE  INJECTION  FACILITY.   -  -  ^ha   operational 
system  employs  two  injection  facilities,  one  in  California 
and  one  in  Minnesota.   Each  injection  facility  contains 
a  60  foot  parabolic  directional  antenna  in  conjunction 
with  a  10  kilowatt  transmitter.   A  small  computer  similar 
to  that  used  in  the  A.N/BRN-3  submarine  navigation  equipment 
controls  the  operation.   The  "injection  mess  ge  from  the 
operations/computer  center  is  received  at  the  injection 
facility  and  stored.    Len  the  satellite  is  in  range  of  the 
injection  facility  the  message  proper  containing  bits  of 
information  in  combinations  of  binary  on^^    and  zeros  is 
converted  to  waveform  to  modulate  the  carrier  wave  of  the 
transmitter.   The  messape  is  transmitted  into  the  satelli^. 
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in  15  seconds.   Since  the  satellite  and  the  earth  are  in 
'notion  during  transmission,  the  distance  from  the  transmitter 
to  the  satellite  varies  constantly.   Therefore,  the  time 
reauired  for  the  signal  to  reach  the  satellite  must  be 
competed  in  advance  and  the  bit  rate  must  be  adjusted  to 
ensure  that  the  satellite  receives  them  at  a  3ven  rate. 
The  bit  rate  is  changed  32  times  during  the  15  second 
transmission.   The  injection  facility  also  reads  the  message 
as  retransmitted  from  the  satellite  to  determine  bhe  success 
of  tne  injection  attempt.   If  necessary,  the  injection 
facility  recycles  to  do  bhe  whole  operation  over  again, 
complete  c^cle  for  each  injection  requires  two  minutes. 
The  equipments  of  the  injection  facility  consist  of 
commercially  available  components  wherever  possible.   In 
fact  the  transmitter  is  a  modified  commercial  television 
transmitter. 

4.6   NAVIGATION  EQUIPMENT.   -  -   Two  types  of  shipboard 
navigation  equipments  have  been  developed,  one  for  submarines 
and  the  other  Tor  surface  ships.   The  AN/BRi\T--3  tvas  developed 
by  Westinghouse  and  is  designed  specifically  for  submarine 
use.   It  incorporates  many  features  adapted  to  submarine 
environment,  and  certain  pieces  of  interface  equipment  to 
permit  operation  in  conjunction  with  other  navigational 
equipments  presently  used  aboard  Fleet  Ballistic  i  issile 

submarines.   The  AN/BRN-3  is  completely  automatic  in 
operation.   It  computes  it's  own  alerts  of  wh  n  navigational 
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satellites  will  bo  within  radio  range  of  the  submarine 
position.   It  then  computes  the  satellite  position  for 
specific  intervals  during  a  pass  based  upon  knowledge  it 
already  has  of  the  satellite  orbit.   Using  an  assumed 
submarine  position  it  precomputes  the  expected  doppler 
curve.   When  the  actual  doppler  curve  is  received  the 
AN/BRN-3  adjusts  and  readjusts  the  assumed  position  until 
the  closest  possible  matching  of  the  two  curves  is  achieved. 
Since  the  assumed  position  is  supplied  by  an  Inert ial  or 
some  other  high  quality  positioning  system,  the  curve- 
fitting  process  takes  place  almost  in  real  time.   The 
position  at  which  the  two  curves  fit  most  closely  together 
is  taken  to  be  the  submarine's  oo^ition.   This  fix  is 
transmitted  to  the  other  navigation  equipment  on  board,  as 
necessary.   In  addition,  the  eouipmont  performs  its  own 
readiness  checks,  diagnoses  any  trouble,  and  {  "ints  out  in 
english  what  part  to  replace. 

The  equipment  developed  for  surface  ships  is  designs. id 

the  AN/SRN-9.   It  is  far  less  complex  than  the  submarine 
equipment  since  it  does  not  have  the  environmental  and 

interface  problems  associated  with  undersea  ooerat:ions. 
The  AN/SRN-9  receives  the  signal  from  tne  satellite  and 
automa tically  computer  the  ship  position.   It  uses  the 
orbital  parameters  provided  by  the  satellite  to  compute 
the  part  of  the  orbit  corresponding  to  the  dooplor  curve 
received.   From  this  it  computes  the  satellite  coordinates 
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at  the  observation  time  points.   ViTith  this  information, 

plus  the  integrated  floppier  c^cle  count  for  at  least  three 

2  minute  time  intervals,  it  computes  the  corrections  dX  , 

dY  ,  and  dZ   to  the  assumed  position  coordinates.   The  new 
n*       n  * 

position  coordinates  are  then  used  to  calculate  the  navigator 
latitude  and  longitude. 
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5.   SYSTEM  ACCURACY 

5.1   GENERAL.   -  -   In  the  course  of  the  development, 
of  the  Navy  Navigation  Satellite  system,  individual  er^or 
components  were  analyzed  in  order  to  determine  and  improve 
those  components  which  most  affect  the  system.   Theoretical 
values  based  on  assumed  idealized  conditions  were  determined 
for  most  components.   These  values  can  be  expected  to  vary 
under  operating  conditions.   Specific  values  pertaining 
to  the  accuracy  of  fix  determinations  based  on  the  use  of 
the  prototype  and  operational  systems  are  not  available 
due  to  security  restrictions.   However,  several  fleet 
units  have  made  extensive  use  of  the  operational  system 
and  fix  accuracies  of  i  0.1  mile  or  better  have  been  made 
public . 

From  the  navigator's  point  of  view  the  primary  error 
components  of  the  operational  system  are: 

1.  Ionospheric  refraction  error. 

2.  Frequency  error. 

3.  Time  error. 

4.  Error  in  the  determination  of  the  velocity 
of  eleccromagnetic  propagation. 

5.  Navigator  instrument  error. 

6.  Navigator  velocity  er^or. 

7.  Satellite  coordinate  error. 

A  short  discussion  of  these  error  component^  is  given  in 
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the  following  sections.   It  should  bo  noted  that  flue  to 
the  nature  of  the  system,  the  satellite  coordinate  errc1 
can  he  broken  down  into  individual  er*ror  components  similar 
to  those  encountered  in  fix  determination.   The  principal 
difference  is  in  the  high  order  gravity  coefficients 
required  for  orhit  determination  and  prediction. 

5.2   IONOSPHERIC  REFRACTION.   -  -   In  Section  3.2  it 
was  shown  that  a  first  order  correction  for  ionospheric 
refraction  can  be  determined  by  the  transmission  of  two 
coherent  frequencies  from  the  satellite.   The  resulting 
doppler  shift  still  contains  higher  order  effects  of  the 
series  (l/f  ).   In  order  to  determine  whether  contributions 
higher  than  the  first  order  are  significant,  the  theoretical 
effect  of  refraction  on  the  doppler  shift  was  considered 
in  Reference  9.   Table  5.1  summarizes  the  results  pertaining 
to  upper  bounds  on  the  various  contributions  to  the  refracted 
doppler  shift  for  various  frequencies  within  the  transmission 
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range.   Except  for  the  vacuum  term,  typical  values  are 
probably  a  factor  of  three  to  ten  smaller.   The  first  order 
term  is  proportional  to  the  time  rate  of  change  of  the 
total  number  of  electrons  in  a  unit  cross  section  tube 
extending  from  the  observation  point  to  the  satellite.   This 
is  the  term  that  is  eliminated  when  two  frequency  doppler 
data  is  used.   The  second  order  term  considers  Faraday 
rotation  and  its  sign  denends  on   whether  the  antenna  is  left 
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or  right  circulary  polarized.   The  third  order  correction 
is  composed  of  four  terms.   The  first  term  arises  from 
the  fourth  order  term  in  the  expansion  of  the  refractive 
index  series.   The  second  term  reoresents  the  difference 
between  the  optical  path  and  the  geometric  path.   The  last 
two  terms  represent  modifications  to  the  solution  o"°  the 
wave  equation  using  the  geometrical  optics  approximation. 
These  terms  depend  upon  various  derivatives  of  the  electron 
density.   Unless  the  satellite  frequencies  are  below  VHP 
the  two  terms  not  obtain;ible  from  the  geometrical  optics 
approximation  are  negligible. 

Table  5.2  presents  the  estimated  maxima  of  the  second 
and  third  order  contributions  to  the  vacuum  doppler  shift 
when  dual  frequency  data  are  combined  to  eliminate  the 
first  order  refraction  term.   Three  typical  frequency 
combinations  are  given  and  the  table  entries  have  been 
converted  from  frequency  shift  to  errors  in  the  r---nge  rate 
in  meters  per  second. 

Examination  of  the  tables  indicates  that  if  the  lower 
of  the  two  frequencies  is  not  less  than  100  mc/s  and  the 
ionosphere  is  not  extremely  disturbed,  the  second  and  third 
order  refraction  contributions  should  be  negligible. 

5.3  FREQU  1NCY  AND  TIMS  ERRORS.  -  -  The  navigator 
must  have  an  accurate  clock  and  a  reference  frequency  in 
order  to  determine  his  position  from  doppler  observations. 
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In  addition  he  muft  know  the  freouency  transmitted  by  the 
satellite  and  the  satellite  position  as  a  function  of  time. 
In  Chapter  3  it  was  shown  that  to  the  requirDd  accuracy  the 
frequency  of  the  reference  oscillator  is  assumed  unknown 
and  the  position  determination  computations  were  arrange^ 
so  as  to  eliminate  the  need  for  this  vnlue.   The  navigator 
is  able  to  recover  time  from  each  satellite  observation  to 
an  accuracy  of  better  than  one  millisecond.   This  time  is 
maintained  by  a  clock  derived  from  the  reference  oscillator 
which  has  a  drift  rate  of  less  than  1  part  in  10    per  day. 
Therefore,  the  time  and  f^oauency  errors  in  the  Navy 
Navigation  Satellite  System  both  derive  from  the  drift  of  the 
satellite  oscillator.   A  drifting  oscillator  gives  an  along 
track  navigation  error  because  the  frequency  is  in  error 
and  causes  a  time  error  in  the  clock  that'  it  is  driving. 
This  time  error   also  gives  an  along  track  error.   Navigation 
errors  of  about  1  kilometer  result  from  an  error  of  1  pr.rt, 
in  10   in  the  satellite  frequency.   However,  satellite 
frequency  can  be  determined  with  an  accuracy  of  at  least 

1  part  in  10    and  predicted  with  an  accuracy  of  at  least 

2  parts  in  10    per  day.   Since  the  satellite  oscillator 
can  be  tuned  by  command  from  an  injection  station  twice 


per    day,    it's    frenuency    can  be   kept    to    a    predetermined 

Go] 


valije  within  an  accuracy  of  a  least  2  parts  in  10   .   This 


causes  a  navigation  error  of  approximately  20  meters 
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5.4  V«LOCITY  OF  EL33CTROMAG  ETIO  PROPAGATION.  -  -   The 
value  for  the  velocity  of  electromagnetic  energy  propagation 
i «?  taken  to  he  a  constant  in  calculations  for  ;">os*tion 
determination  using  the  ?T,^vy  Navigation  Satellite  ^y^tem. 
The  value  is  probably  known  to  an  accuracy  of  1  part  in 
3,000,000.   This  would  oermit  position  determination  to  on 

GO 

accuracy  of  about  two  meters.  "   Future  improvements  in  this 
constant  can  be  accomodated  by  modification  of  the  system 
based  on  the  current  accepted  value. 

5.5  NAVIGATION  INSTRUMENT  ERROR.  -  -   In  order  to  test 
the  contribution  that  the  navigation  equipment  rakes  to  the 
navigation  error,  observations  were  made  on  experimental 
satellites  transmitting  on  four  frequencies.   Two  navigation 
instruments  at  the  same  site  made  simultaneous  observations 
on  the  satellite.   One  instrument  used  two  of  the  satellite 
frequencies  and  the  other  used  the  remaining  two  frequencies. 
The  difference  between  the  navigated  positions  is  then  the 
result  oT  two  independent  measurements  of  position.   The 
result  of  one  set  of  tests  obtained  with  the  Transit  2A 
satellite  gave  a  mean  difference  in  position  of  50  meters 

E3 

m  one  coordinate  and  30  meters  in  the  other.    It  is 
expected  that  instrumentation  hap  improved  since  this  time. 

5.6  NAVIGATOR  VELOCITY  ERROR.  -  -   Since  the  navigator 
must  use  a  measured  value  of  hi^  velocity  to  determine  his 
position  at  the  times  of  the  satellite  observations,  an 
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error  in  this  m^afTement  will  cause  on  err-or  in  the  fix.   A 
ship  acted  on  by  unknown  current,  wind  ope  drift,  and  pes 
condition  effects  amounting  to  a  one  knot  velocity  error  in 
the  northerly  direction  has  an  estimated  position  error  of 
0.3  nautical  miles.  "   A  one  knot  easterly  velocity  error 
causes  a  negligible  positional  error  at  mid-latitudes  since 
the  satellite  is  in  a  polar  orbit.   It  was  determined  to  at 

least  three  olace  accuracy  that  the  position  errors  varied 

ft) 

linearly  with  the  velocity  errors.    This  does  not  appear 

to  be  a  severe  limitation  for  surface  ships  and  submarines 
where  accurate  pit-lop  or  inertial  system  velocities  are 
normally  available. 

5.7   SATELLITE  COORDINATE  ERROR.  -  -   Apart  from  the 
oreceedinp  error  components  which  are  normally  negligible, 
the  only  remaining  source  of  system  er^or  is  in  the  satellite 
coordinates  transmitted  to  the  user.   Any  error  in  these 
coordinate^  is  reflected  almost  exactly  one  for  one  in  the 
accuracy  of  the  navigation  fix.   Tests  on  the  size  of  the 
errors  in  the  predicted  orbits  of  the  Navy  Navigation 
Satellite  System  were  conducted  by  computing  and  predicting 
the  orbits  at  four  tracking  stations.   These  predictions 
were  made  using  geodetic  parameters  in  use  at  the  Naval 
Weapons  Laboratory  during  four  periods  of  time  through  1964. 
The  predicted  orbit  was  compared  with  an  orbit  ba^ed  upon 
independent  doopler  observation?  made  by  world-wide  tracking 
stations  during  the  prediction  interval  and  using  the  latest 
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available  gravity  coefficients.   The  resulting  twelve 
hour  prediction  error?  for  random  set?  of  observations 
made  on  polar  orbiting  satellites  show  o  decrease  in 
size  from  over  1000  meters  in  1960  to  the  value  of  50 
meters  in  1965. 
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5.8  PROPOSED  ERROR  ANALYSIS.  -  -  A  high  speed  computer 
may  be  easily  programed  to  analyze  the  effect  of  individual 
and/or  collective  error  component?  on  a  given  satellite 
position  determination.   The  basic  inputs  of  this  program 
would  consist  of  the  following  parameters; 

1'      V  V  Zn 

2.  iAx,  iAy,  i/z 

3.  X^,  Y^,  Z£ 

4.  fr 

5'   fs 

6 .  c 

7.  AT 

The  satellite  coordinates  as  a  function  of  time  in  an 
earth  centered  rectangular  coordinate  system  are  available 
from  a  Geodetic  Science  Computer  Library  Program  at  the 
Ohio  State  TTniversity  or  from  some  other  similar  program. 
TTsing  satellite  coordinates  bracketing  three  or  more  time 
intervals  and  choosing  the  above  parameters  which  are  assumed 
errorless,  we  can  determine  the  expected  integrated  cycle 
count  from  (3.7).   Assuming  errorless  instrumental  and 
environmental  conditions  the  calculated  cycle  counts  and 
the  above  oarameters  may  be  considered  as  an  ideal  errorless 
position  determination  based  on  satellite  observations. 
This  can  be  verified  by  substitution  of  these  values  into 
the  equations  developed  in  Chapter  3.   The  position 
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coordinates  produced  by  calculation  will  be  equal  to  the 

initial  navigator  coordinates  since  the  corrections  dX  , 

dY  ,  and  d/,f   will  be  zero  in  our  assumed  errorless 
n7         n 

condition.   We  have  thus  created  an  ideal  errorless 
position  determination  which  can  be  used  as  a  standard 
in  an  error  analysis.   By  random  and/or  systematic 
alteration  of  the  assumed  errorless  input  parameters  we 
can  determine  their  individual  and/or  combined  effect  on 
the  final  position  determination.   A  study  could  also  be 
made  to  determine  the  effect  of  decreasing  the  time  interval 
and  using  more  observations  in  a  least  souares  solution. 
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6 .   SUMMARY 

Tho  adaption  of  space  technology  provides  a  new  aporoach 
to  the  earth-bound  problem  of  marine  navigation.   The  Navy 
Navigation  Satellite  System  is  based  on  the  floppier  shift 
in  a  radio  frequency  transmitted  from  a  near-earth  satellite. 
By  exact  measurement  of  the  doppler  shift  it  is  possible  to 
determine  the  location  of  the  receiver  on  the  surface  of  the 
assumed  figure  of  the  earth.   The  accuracy  obtained  by  using 
this  technique  is  possible  because  the  satellite  is  very 
precisely  restricted  to  an  orbit  defined  by  physical  laws 
of  motion  and  the  quantities  observed,  time  and  freauency, 
can  be  measured  to  an  accuracy  of  one  part  in  10   .   Of  all 
the  possible  satellite  paths  permitted  there  is  only  one 
which  results  in  a  particular  curve  of  doppler  shift. 

Two  basic  problems  are  inherent  in  this  system.   One 
is  that  the  iono^ohere  bends  radio  waves  and  thus  gives  a 
false  position  of  the  satellite.   This  has  been  overcome 
by  the  use  of  two  satellite  transmission  frequencies  which 
allow  a  first  order  correction  to  be  computed.   The  second 
basic  problem  is  the  determination  of  the  gravity  field  of 
the  earth.   This  affects  the  determination  and  prediction 
of  satellite  orbits.   This  problem  has  been  reduced  to  an 
acceptable  level  by  using  satellite  orbital  data  to  calculate 
new  gravity  coefficients. 

The  use  of  a  navigational  satellite  system  based  on 
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the  dopoler  shift  of  VHF  radio  waves  appears  to  o^ovide 
the  best  oresent  day  solution  to  the  military  navigation 
requirements.   World-wide  coverage  i^  obtained  by  use  of 
polar  orbits.   The  system  is  inherently  all  weather  and 

immune  to  interference  since  line-of-sipht  radio  frequencies 
are  used  for  both  tracking  and  navigation.   ^ince  the 
satellites  transmit  their  ephemeridas  the  navigator  need 
neither  interrogate  the  satellite  nor  receive  orbital  and 
time  information  by  other  communication  links.   The  system 
is  also  reliable  because  a  temporary  interruption  of 
reception  during  a  pass  does  not  preclude  the  use  of 
data  that  are  received  to  obtain  an  accurate  navigational 
fix.   It  is  assumed  that  the  system  meets  the  military 
accuracy  reauirements . 
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